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In the covariant light-front quark model, we investigate the Be hc,Xc0.i,2 form factors. The 
form factors are evaluated in space-like kinematic region and are recasted to the physical region by 
I adopting the exponential parametrization. We also study the semileptonic Be decays and find that 

' branching fractions for the Be {he, Xc0,i,2)li^{l ~ e,^) decays have the order 10"'^ while branching 

' fractions for Be — > {he, Xc0.i,2)tVt are suppressed by one order. These predictions will be tested on 

j-H . the forthcoming hadron colliders. 

^ : 

'"J . PACS numbers: 13.20.He, 12.39.Ki 

^ : 

(N . 

^ ^ I. INTRODUCTION 

^ ■ 

I ' B meson weak decays provide a golden place to extract magnitudes and phases of the Cabibbo-Kobayashi-Maskawa 
Oh' 

^ ' (CKM) matrix elements, which can test the origins for CP violation in and beyond the standard model (SM). 
i ' Semileptonic and nonlcptonic B meson decays have received extensive interests and achieved many great successes. 
[ Experimentally, the two B factories have accumulated more than 10^ B — B events; measurements are becoming more 
^ ' and more precise. On the theoretical side, apart from contributions proportional to the form factors, the so-called 
' nonfactorizable diagrams and some other radiative corrections are also taken into account. All of these are making B 
] physics suitable for search of new particles and new phenomena (see Ref. [ll for a review) . 
• I Compared with B decays, B^ meson decays have received much less experimental considerations. The mass of 
' a BcBc pair has exceeded the threshold of T(4S') thus the Be meson can not be generated on the two B factories. 
But Be meson decays have a promising prospect on the forthcoming hadron colliders. The Large Hadron Collider 
^ ^ (LHC) experiment, which is scheduled to run in the very near future, will produce plenty of B^. events. The LHCb 
collaboration has the desire to perform a comprehensive investigation on B^ meson decays. With more and more data 
5-H ' accumulated in the future, the study on Be mesons will be of great interests: (1) Be contains two different heavy 
flavors, the spectroscopy may be different with the light meson or the meson with only one heavy quark. It serves as 
a different laboratory to study the strong interactions. (2) Be meson can weakly decay via the b ^ q transition like 
lighter Bu,d,s mesons, but the dynamics is dramatically different. (3) Moreover, the charm quark can also decay via 
weak interactions, where the b quark acts as a spectator. The CKM matrix element \Ves\ '--^ 1 is much larger than 
\Veb\ ~ 0.04 in h quark decays. Decays of the charm quark contribute much more to the decay width of the Be meson. 
Although the phase space in c — )■ d, s decays is much smaller than that in 6 ^ c transitions, the former ones provide 
about 70% contributions to the decay width of Be- This results in a larger decay width and a much smaller lifetime 
than the B meson: tb^ < ^tb- (4) The two heavy h and c quarks can annihilate and they provide new kinds of weak 
decays with sizable partial decay widths. The pureJe_£tonic or radiative leptonic decay can be used to extract the Be 
decay constant and the CKM matrix element Veb 

Semileptonic Be decays are simpler than nonleptonic Be decays: the leptonic part can be perturbatively evaluated 
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leaving only hadronic form factors unknown. In two-body nonleptonic decays, most channels are also dominated 
by the Be transition form factors. Thus the study of Be transition form factors is essentially required. In the present 
work, we will use the light-front quark model to analyze the form factors Be decays into p-wave chamonia. This can be 
viewed as a continuation of our previous work P| . The light front QCD approach has some unique features which are 
particularly suitable to describe a hadronic bound state [5[- Based on this approach, a light-front quark model with 



many advantages is developed 0, 0, S, [s , llfll • This model can provide a relativistic treatment of the movement of the 
hadron. Is also gives a fully treatment of the hadron spin by using the so-called Melosh rotation. The light front wave 
functions are expressed in terms of their fundamental quark and gluon degrees of freedom. They are independent of 
the hadron momentum and thus are explicitly Lorcntz invariant. In the covariant light- front quark model, the spurious 
contribution, which is dependent on the orientation of the light-front, becomes irrelevant in physical observbles and 
that makes the light-front quark model more selfconsistent. This model has been successfully extended to investigate 
the decay constants and form factors of the s-wave and p-wave mesons [ll|, Q, 13 1. 



Our paper is organized as follows. The formalism of the covariant light-front quark model and numerical results for 
form factors are presented in the next section. The decay rates of semi-leptonic B^ decays are discussed in Section Hill 
Our conclusions are given in Section HVl 



II. FORM FACTORS IN THE COVARIANT LIGHT-FRONT QUARK MODEL 

Be ^ S, A (S, A denotes a scalar meson or an axial- vector meson, respectively) form factors are defined by 

/ 2 2 \ 2 2 

{SiP")\A,\BeiP')) = i^P, ""^ g^j Ff^'iq^) + "^' q.F.^^'iq^ (1) 

{AiP",e"*)\WBeiP')) = e^.o/3£"*''F"/A^^^(<z2), (2) 

niB, - mA 

e"* ■ P 

m-B, - niA 

e"^P 



{A{P",e"*)\V,\Be{P')) - -^{{mB^-mA)s';,*Vf^^{q') P.V.'^^^iq') 

-2mA '—^ [Vf^-^iq') ~ ^o'^^^l?^)] }, (3) 



where P = P' + P" , q = P' — P" and the convention eoi23 = 1 is adopted. To smear the singularity at = 0, the 
relation Vf'^iO) = Vff"^(0) is reqmred, and 
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Form factors of Be decays into a tensor meson are defined by 

{T{P",e"*)\V^\Be{P')) = h{q^)e^,e.t5e"*''^PxP''q^ 

{T{P'\e"*)\A,\Be{P')) = ~i{k{q^)ellP'^+e%P-P^[P,h+{q^)+q,h^{q^)]], (5) 

where the polarization tensor, which satisfies tp^uP"^ = 0, is symmetric and traceless. The spin-2 polarization tensors 
can be constructed using spin-1 polarization vectors: 

6V(P,±2) = e^(±)e.(±), eV(p,±l) = -i=[e^(±)e,(0) + e,(±)e^(0)], 

e^-b^O) = -L[e^(+)e,(-) + ^e^{0)e,{Q), (6) 
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where e is the polarization vector for a vector meson. If the recoiUng meson is moving on the plus direction of the z 
axis, their explicit structures are chosen as 

£^(0) = — (Ipt 1,0,0, i;T), 



1 



e^(±) = _(0,Tl,-i,0), (7) 



where the Et and is the energy and the magnitude of the momentum of the tensor meson in the rest frame, 
respectively. rriT denotes the tensor meson's mass. In B^, meson decays, it is useful to define a new polarization vector 
ct for the tensor meson 

eTfj.{h) =€fj,,y{p,h)pB^ (8) 

which satisfies 

eT^(±2)=0, eT^{±l) = ^eiO)-pBM±)^ ^T^^iO) ^ ^e(0) ■ psM^)- (9) 

The contraction is evaluated as e(0) -pB^ = mBc\'PT\/mT- We can see that the new polarization vector plays a similar 
role with the polarization vector for a vector meson, regardless of the nontrivial factors ^ or In analogy with 
the Bc^ A transition, one can define the following form factors for convenience: 

A''^^ = -{mB^-mT)h{q^), Vf^^ = V^^^'' = {mB^ - mT)h+{q% 

^ ^V-W) - '^vr-ii) - ^^-(,^ (10) 

where these form factors A^"'^ , ^^1,2 have nonzero mass dimensions. 

We will work in the g'+ = frame and employ the light-front decomposition of the momentum P' = (P'~, P'+, P^), 
where P'^ = P'^ ± P'^, so that P'^ = P'+P'~ — P'^. The incoming (outgoing) meson have the momentum of 
P' = p'l +P2 ( P" = p'i +P2) and the mass of M' (M"). The quark and antiquark inside the incoming (outgoing) 
meson have the mass m'/"^ and 77x2 and the momenta are denoted as p'/"^ and P2 respectively. These momenta can 
be expressed in terms of the internal variables {xi,p'j_) as: 

p'+ = a;i,2P'+, p'l, 2± = a;i,2Pi ± p'^, (11) 



with Xi + a;2 = 1. Using these internal variables, one can define some useful quantities for the incoming meson: 

"'2^r< ^ pf_ 

Xi X2 



= (e;+e2)2 = ^i±^ + ^i±^, K = JMi^-{m\-m2Y 



here can be interpreted as the energy of the quark or the antiquark and M'q can be viewed kinematic invariant 
mass of the meson system. The definition of the internal quantities for the outgoing meson is similar. To calculate 
the amplitude for the transition form factor, we require the following Feynman rules for the meson-quark-antiquark 
vertices (iF^): 

iV'p = H'pjr,, 

iT's = H's, 
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FIG. 1: Feynman diagram for Be S,A,T decay amplitudes. The X in the diagram denotes the V, A transition vertex while 
the meson-quark-antiquark vertices are given in the text. 

iV'rp = i^H^bf,- r^ip'l-P2)^]{p'i-p2)u- (13) 

For the outgoing meson, one should use i(7or2/7o) for the corresponding vertices. 

In the conventional light-front quark model, the constituent quarks are required to be on mass shell and the physical 
quantities can be extracted from the plus component of the corresponding current matrix elements. However, this 
framework suffers the problem of non-covariance because of the missing zero-mode contributions. In order to solve 
this problem, Jaus proposed the covariant light-front approach which permits a systematical way to deal with the 
zero-mode contributions [lo| . Physical quantities such as decay constants and form factors can be calculated in terms 
of Feynman momentum loop integrals which are manifestly covariant. The lowest order contribution to a form factor 
is depicted in Fig. [1] The X in the diagram denotes the V, A transition vertex. In Be to p-wave chomonia decays, 
Pi{Pi) is the momentum of the bottom (charm) quark, while p2 is the momentum of the antiquark. It is also similar 
for the notation of the quark masses. As an example, we will derive the Be S transition amplitude 

^ " ' (2^)4 y "^^^N'^N'lm ^ ' ^^^^ 
where N'^"^ = p'/""*^ — to'/"''^, N2= p\ — Tn\. The function S^"^ is derived using the Lorentz contraction 

5f ^ - Tr [(// + m'/)7A.75(/i + i>2 + ^2)] 

= 2p'i^[M'2 + - g2 _ 2iV2 - (m'l - m^f - {rn![ + m^f + (m'l + m'/)^] 
+q^Xq^ - 2Af'2 + iV{ - N'i + 2N2 + 2(m'i - m^f - (m'l -f m![f\ 

+P^{q^ N'i - (m'l + m'/)']. (15) 

In practice, we use the light-front decomposition of the loop momentum and perform the integration over the 
minus component using the contour method. If the covariant vertex functions are not singular when performing the 
integration, the transition amplitude will pick up the singularities in the antiquark propagator. The integration then 
leads to 
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tt' tt" qB^S 



where 



112 I 112 

112 _ Pl_ +^1 



Xi X2 

with p'^ — p'j_ — X2q±- The exphcit forms of h'j^j and ui^^ used in this work are given by 



XiX2 1 , 



Ml,' 



7^P 



(16) 

(17) 

(18) 
(19) 



-, WiA — A 



X\X2 1 , 

Nc ViK^^" 



m\ — 1712 



(20) 



where ip' and is the hght-front wave function for s-wave and p-pave mesons, respectively. After this integration, 
the conventional light-front model is recovered but manifestly the covariance is lost as it receives additional spurious 
contributions proportional to the lightlike four vector uj = (0, 2,0j^). The undesired spurious contributions can be 
eliminated by the inclusion of the zero mode contribution which amounts to performing t he p ~ integration in a proper 
way in this a 
Appendix A. 



way in this approach. The specific rules under this p integration are derived in Ref. 
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and are displayed in 



Using Eqs. (fT51) - (pn)) and taking the advantage of the rules in Ref. [lOljlllii the Be ^ S form factors are straight- 
forwardly given by 

Nc 



h' h" 

dx2Cpp'^- 



Xi{M'^' + M'^') + X2q' 



167r3j - '-X2N[N'{ 

-X2{m[ + m'lY — xi{m'i — 7712)^ — Xi{m'l + 7712)^ 



Nc 

q ■ P 167r 



3 ' X2N[N'{ 



X1X2M' — p1 — m[m2 — {ni'l + ni2){x2m[ + Xim2) 



q^ \ q^ J 9" 

-{x2 - a;i)M'2 -f 2xiM'^ - 2{m[ - 7772)(mi - 777") 



M" 



X2{q' + q-P) 



(21) 



Similarly, one can derive the Be A,T form factors and we refer to Appendix B for tedious expressions of these form 
factors. 

The light front wave function (p' can be obtained by solving the relativistic Schrodinger equation with a phe- 
nomcnological potential. But in fact except for some limited cases, the exact solution is not obtainable. In practice, a 
phenomenological wave function to describe the hadronic structure is preferred. In this work, we will use the simple 
Gaussian-type wave function which has been extensively examined in the literature [lli ll^ . Iisl 

3/4 



, , , , . ri~ , dp'^ e[e2 



(22) 
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The parameters (3's, which describe the momentum distribution, are usually fixed by mesons' decay constants whose 
analytic expressions are also given in The decay constant for the Be meson is employed by 



/i3^ (400±40)MeV, (23) 
which gives Pb^ = (0-89jlQ Qy4) GeV. This value is a bit smaller than results provided by Lattice QCD method 

Jb^ = (489 ± 4 ± 3)MeV. (24) 



The other inputSjincluding masses (in units of GeV) of the constituent quarks and hadrons, Vet and the lifetime of 
Be, are used as [l5| 



1.4, rub = 4.8, niB^ 



6.286, 



Vcb- 



m 



= 3.52528, Tn^^„ = 3.41476, 
41.2 X 10"^ TB, = (0.46 ± Om)ps. 



3.51066, 



''Xc2 



3.5562 



(25) 



13l |. The shape parameter for Xci is 



The constituent quark masses are close to those used in the literature 
used as: (3^^-^ ~ (G.7±0.1) GeV which corresponds to \fx^i \ = (340ti[J?) MeV. For the other shape parameters, we will 
assume the same values and introduce a relatively large uncertainty to compensate the different Lorentz structures: 
/3xco = /3xc2 = = (0.7 ±0.1) GeV. 

Unlike the light quark, the heavy bottom and charm quark have large masses. In the heavy quark limit m^^b oo, 
the Be and charmonium systems obey the heavy quark symmetry which is helpful to simplify the dynamics in transition 
form factors and decay amplitudes. In particular, the large momentum of the heavy quark can be projected out and 
the remanent momentum is of the order of the hadronic scale. For Be and charmonia, the two constituents are both 
heavy and move non-relativistically. After projecting the large mass scale, the dynamic scale is of the order nieV 
and niev'^, where v is the relative velocity of the quark-anti-quark pair. Then physical quantities can be expanded 
in terms of l/nie, l/rrib in the effective theory. In the present analysis, the expansion in l/nie and l/nib is not used 
and physical quantities contain a tower of contributions with different orders. The Be and the charmonia are directly 
made of two heavy quarks and the dynamics is reflected by the light-front wave functions. Despite of the different 
treatments, the leading power behavior should be the same. In Be and charmonia, the transverse momentum is of 
the order of meV. Inferred from Eq 
In nonrclativistic QCD (NRQCD) 



22|) . the parameter /?' in the light- front wave function is of the order of meV. 
16|, the kinematic energy has the order of iriev'^. If it can be identified as the 



hadronic scale , the /?' is of the order of y/mc^QCB, which is expected to be larger than the shape parameter in the 
light meson system. This feature is also confirmed by the numerical result: the shape parameters /3's for charmonium 
and Be meson are larger than those for light mesons such as = 0.3102 GeV 11 1. 



In the NRQCD framework, the light-cone distribution amplitudes of the s-wave charmonia have been comprehen- 
sively investigated in Refs. 



19 



20| . Recently, the analysis has been generalized to the p-wave charmonia [2l[ . In 



their notation, the distribution is described by the matrix elements of the nonlocal operators, while in the light-front 
quark model, we use the coupling vertex. Moreover, the distribution amplitudes in these two frameworks are also 
different. In Ref. 2]|, the leading twist light-cone distribution amplitude is expanded into Gegenbauer polynomials 



and the Gegenbauer moments are studied in the QCD sum rules. The authors in Ref. 2l| also propose the following 
model 



$(^,/i^m,) = c(/3p)(l-e')^exp 



f3f 



1-e 



(26) 



7 



TABLE I: Results for the Be XcO,ci,c2,^c form factors and fitted parameters ci and C2. The first type of uncertainties are 
from the shape parameters of the p-wave charmonia and the second ones are from the Be meson decay constants. 



F 


F(0) 




Cl 


C2 


rpBcXcO 


'^•^'-0.06-0.01 
n ^7+0.03+0.00 
-0.05-0.01 


r| 70 + 0.01 + 0.02 
' 'J-0.05-0.03 

n /in+0-00+0'00 

'^•^'^-0.02-0.01 


1 A A-? 1 A 1 A 

9 riq+0.2 ^ +0. 10 
Z.UO_o 24-0.10 
n /|r+0.33+0.00 
'-'•^^-0. 29-0. 00 


n /1 0+0. 05+0. 01 

'^•^■J-O.06-0.01 
1 01 +0.12+0.09 
J^-'J^-O.ll-O.lO 


^0 

yB.h^ 


r. r,7+0.00+0.01 
'J-"' -0.01-0.01 
r, f;^+0.10+0.02 
'-'•"^-0.00-0.02 
n cn+0.03+0.04 
U-OU.Q 06-0.05 
r, 09+0.05+0.03 
"•■J^-0.04-0.03 


r, 1 1 +0.00 + 0.01 

'"'■^^-0.01-0.01 

n Q/l+O.03+0.03 
^■^^-0. 08-0. 04 
n c;+o.Oi+u.04 
^■">J-0.05-0.05 
n cc+0.06+0.06 
U-'JO_o.03-0.05 


9 09+0.30+0.11 
^"^^-0.26-0.11 
1 Q9+0.35+0.04 
^■='^-0.32-0.05 

1 A OA 1 A 1 1 

1 c:/l+0.2y+o.il 
-•-■^^-0.24-0.11 
9 f;q+0.28+0.09 
Z.UO_o 26-0.10 


r, 4Q+0.07+0.01 

'^•^■'-0. 07-0. 02 
r, qQ+0.08+0.00 
"•■J ="-0.10-0. 01 

1 A A(^ 1 A A1 

n 9/1+0. U6+0. 01 
"•^^-0.06-0.02 

n fiq+o.06+0.01 

"•"•3-0.07-0.01 


j^B^Xal 
T/BcXcl 
yBaXal 
yBaXal 


n qe+0.02+0.01 
U.OU_Q o4_o 02 

1 q + 0.01 + 0.00 
'-'•J^'J-0.01-0.01 
^ c^^+n 00+0 02 
U.OO_Q Q2_o.o3 

n 1 c+o.oi+o.oi 

'-'•^^-0.01-0.01 


n c/1+0.01+0.03 

'J-'-'*-0.03-0.04 
r, 90 + 0.01 + 0.00 
^■^"J-0.02-0.00 
r\ 70+0. 05+0. 03 
''J-0.06-0.04 

n 1Q+0.00+0.02 

^■^^-0.01-0.02 


1 QS+0-26+0.11 
l.»0_o 23-0.11 

9 QQ+0.19+0.32 
^■^="-0.13-0.29 
^ r T +0 40+0 10 
>J-jr_o 35_o.io 

1 99+0.18+0.20 
^•^^-0.16-0.20 


n /,q+o.o4+o.oi 

"•l'J-0.06-0.02 
n 090+0.32+0.16 

U.UZO_J^ Q4_Q 3g 

OQ+0. 19+0.01 
— i.rf5_o i7_o.o2 

n nc+0.26+0.09 
U.UO_o 05-0.11 


j^ScXc2 
^ScXe2 
^ScXe2 


r, r>99+0. 002+0. 001 
"•'-'^^-0.003-0.001 
1 97+0.15+0.05 
J^-^ ' -0.13-0.07 

n m 1 +0.001+0.000 

U.U11_Q 000-0.000 

n non+o-oo+o-oo 

U.UZU_o 00-0.00 


n n'^R+0-002+0.003 

U.UOU_o 004-0.003 

1 70+0.12+0.11 

''J-0.25-0.14 

n ni c+0. 002+0. 000 

U.UiO_o 001-0.000 

n n'?'?+0-002+0.003 

U.UOO„o. 004-0. 003 


9 ro+0.28+0.01 
z.oo_o 25-0.10 
1 f;i +0.21 + 0.12 
^■"^-0.21-0.12 
9 97+0.36+0.07 
-0.41-0.07 
9 /|o+0.26+0.13 
^•^"-0.24-0.11 


n CI +0.06+0.02 
'^•"^-0. 06-0. 01 
r, 9/1+0.01+0.02 
"•^^-0.04-0.02 

n 4f;+0.10+0.01 
"•'±"-0.25-0.01 
n cc+0.05+0.11 
"•iJ"-0.06-0.02 



where c{f3p) is a normalization constant. The parameter ^ is defined as ^ = 2u — 1, where u is the momentum fraction 
of the charm quark. This form is dramatically different with the one in Eq. (|22p used in this framework. The main 
part in Eq. is the exponential term exp (^—jz^^ which corresponds to cxp (^—jpTj^ in the distribution amplitude 
in Eq. . For charmonia, the momentum pf is simplified as 

(x^-x^)^ (27) 
X1X2 

Considering the longitudinal part, we can see that there is an additional factor (xi — 2:2)^ — ^'^ in the one used in this 
framework. It indicates that the distribution amplitude used in this work is sharper at the region around X2 ^ 0.5. 
At last, the parameter /3p is dimensionless and is different with /?'. In the heavy quark limit, u is close to 1/2 and f3p 
is of the order 1: (3p — 3. 4^Qg. Comparing the longitudinal part of the two distribution amplitudes, one can obtain 
the relation between /3' and (3p: f3p ~ 2m^((zi-3:2) ) ^ r^-j^^ typical value of ((a;i — 2:2)^) is of the order ^'^^ which also 
indicates /3p ~ 1. These different distribution amplitudes are expected to induce sizable differences to the resultant 
form factors. They can be discriminated or constrained by the available data of transition form factors in the future. 

Because of the condition = imposed during the course of calculation, form factors can be directly studied only 
at spacelike momentum transfer = ^ which are not relevant for the physical decay processes. It has been 
proposed in llj to parameterize form factors as explicit functions of q^ in the space-like region and then analytically 
extend them to the time-like region. To shed light on the momentum dependence, one needs a specific model to 
parameterize the form factors and we will choose a three-parameter form 

= F(0)exp(cis + C2s2), (28) 

where s — q^ /m\ and F denotes any one of the Be S,A,T form factors. In the fitting procedure, form factors in 
the region q^ = [—1,-15] GeV^ are studied and the fitted results for ci and C2 are collected in table HI We should 
point out that we have adopted a negative decay constant for Xci so that the Be ^ A form factors are both positive. 
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III. SEMILEPTONIC DECAYS 

With the form factors at hand, one can directly perform the analysis of semileptonic B^. decays whose differential 
decay widths are given by 

2 



dr(Bc SIP) 



dr±(_Bc AlP) 
d^ 



1 



38im% 7r3 

Be 



384?Tt| 7r3 



+ (mf + 2g2) 



9 9 



2m/ 



- - q ){mBe 



■^A)V,^^^(q^) 



X(m% m\,q'^) 



■niBe 



■ rriA 



(29) 



(30) 



^A(m|^,m^,q2)G|y- 



384?Tt| TT^ 



X <^ {nq +2q^)\(nil ni\,q^) 



(31) 



where the superscript +(— ) denotes the right-handed (left-handed) states of axial- vector mesons, while the subscript L 
denotes that the axial- vector in the final state is longitudinally polarized, mi is the lepton's mass and X{rri^ , mf, q^) = 
(rrv^ + mf — q^)^ — 4to^ mf with i — S, A. The combined transverse and total differential decay widths are given by: 



dT 
d^ 



dVL 
dq^ 



dVi 



dq^ dq^ dq^ ' dq^ dq^ dq^ 
Expressions for the decay width of Be ^ TlD can be obtained by the decay width of Be 
dTL{B^TlD) _ l\{ml^,m^,q^) dV l{B ^ AID) 



(32) 



Alv decays: 



dq^ 

dT^{B TW) 



Am\ 



dq^ 



2X{m%^,m^,q^) dT^{B ^ 
3 4m^ dq'^ 



AW) 



l(v- 



(33) 



where the form factors Vb,i,2,^ of B^ 

dTL(B^AW) I 



T decays have nontrivial dimensions and thus the two functions 



ventional ^ 



and %^^^^,A^ 

It is compensated by the prefactor 



ABaT\ do uot havc the correct dimensions with the con- 

A(m^ ,m^,q'^) 



which also have nonzero mass dimensions. 



Our results for these semileptonic Be decays are collected in table [TTl Since electrons and muons are very light 
compared with the charm quark, we can safely neglect the masses of these two kinds of leptons. The uncertainties 
are from those in the form factors: the first kind of uncertainties are from the shape parameters of the charmonia and 
the second ones are from the Be decay constant. The third uncertainties are from that in the Be lifetime. Several 
remarks on the results are given in order. First of all, branching fractions for the Be — > {he, Xc0.i,2)li^il = e, ^) decays 
have the order 10""^ while branching fractions of Be — > {he, XcO,i,2)''"f'T are suppressed by one order. In the covariant 
light-front quark model, branching fractions of the Be rjelv and Be J/tplv decays are about one percent Q: 



Bn{Be 



rjeCVe 



(0 57+0 04-1-0. 04-1-0. 10^ 



-0.07-0.04-0.10-' /"i 



Bn{Be Jl^eve 



/I /iQ-l-0.01-l-0.15-l-0.23\q/ 
1^-^^-0.03-0.14-0.23-' 



(34) 



Compared with these decays, the branching ratios of Be — > {he, Xco, 1,2)^2^ are smaller by a factor of 2 — 10. There 
are two main reasons for these differences: the form factors and phase space. For example, if we set the mass of 
Xco equal to that of r\e, the branching ratio of Be ~> Xco''^ becomes 0.40%. The larger form factors of Be — > r\e 
will enhance the branching fraction by a factor of 1.68. Secondly, polarizations of i?c hclv and Be —* Xc\^^ 
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TABLE II: Branching ratios (in units of %) of semileptonic Be — > {hc,XaO,i,2)li^{l = e, and Be — » (hc,XcOA;2)Ti'T decays. 





Be Xcolf 


Be Xcilv 


Be helv 


Be Xc2lf 


This work 


n. oi +0.02+0.01+0.03 
'-'•^-'--0.04-0.01-0.03 


r, 1 <+0. 00+0. 01+0. 02 
'-'•-'-^-0. 01-0. 01-0. 02 


r, n-i +0.05 + 0.01 + 0.05 
'-'•'' -'--0.08-0. 01-0. 05 


r, 1 7+0.04+0.02+0.03 
"•^'-0.06-0.02-0.03 


CCWZ [22] 


0.12 


0.15 


0.18 


0.19 


IKS [23] 


0.17 


0.092 


0.27 


0.17 


IKS [24] 


0.18 


0.098 


0.31 


0.20 


HNV[25] 


0.11 


0.066 


0.17 


0.13 




Be — > XcOTt^r 


Be XcI'TI't 


Be heTUr 


Be Xc2Ti'r 


This work 


"■'-'^^-0.003-0.001-0.004 


n ni p:+o. 000+0. 001+0. 002 

U.UiO_Q QQ^_g QQ2_o.oo2 


"•'-'^^-0.004-0.000-0.003 


n nnQf?+0. 0019+0. 0013+0. 0015 
U.UU»0_oQQ29_o. 0014-0. 0015 


CCWZ [22] 


0.017 


0.024 


0.025 


0.029 


IKS [23] 


0.013 


0.0089 


0.017 


0.0082 


IKS [24] 


0.018 


0.012 


0.027 


0.014 


HNV[25] 


0.013 


0.0072 


0.015 


0.0093 



are dramatically different. As indicated from table HI the form factors Vi and V2 for Be he have different signs. 
Thus the longitudinally polarized decay width receives constructive contributions as we can see in Eq. ([30)1 . The 
form factor A^"'^" is small which suppresses the transversely polarized decay width. Accordingly, a large is 
expected: y~ — ^ hclP decays and ~ 4.7 for Be — > hcTVT- The situation is dramatically different 

for Be Xci decays. The form factors y^"^"^ and vf"^"^ have the same sign, which gives destructive contributions 
to the longitudinally polarized decay width. Form factors A^"^"^ and V^"^"^ are large and thus the minus polarized 
decay width is large. The polarization fraction is reduced and predicted as: ~ 0.24 for Be — > Xcilv decays and 



^ ~ 0.27 for B, 

1 T 



XeiTVr- Compared with results in Refs. [2^, |23|, |2J, |25| for the Be —> {he, Xco, 1.2)11^ decays which 
are collected in table HIl we can see that most of our predictions on the semileptonic Be decays are comparable with 
their predictions. These results will be tested at the ongoing and forthcoming hadron colliders. 



IV. CONCLUSION 



Due to its unique properties, B physics has attracted abundant attentions. Measurements on the CKM matrix 
elements are becoming more and more accurate, which makes the goal to test the CP origins in and beyond SM 
much more practicable. Be meson decays provide another promising place to continue the errand in B meson decays, 
and offer a new window to explore strong interactions. Although the Be meson can not be generated on the two 
B factories, it has a promising prospect on the ongoing and forthcoming hadron colliders. The high statistics of Be 
meson at the forthcoming hadron colliders can compensate for the hadronic pollution and make it suitable for the 
precise determination of many standard model parameters. Because of these interesting features, we have studied the 
Be transition form factors in the covariant light-front quark model, which are relevant for the semileptonic Be decays. 

Branching fractions of Be ~^ {he, XcO,i.2)l'^{l — e,fi) decays have the order 10~^ while branching fractions of 
Be {he, XcO,i.2)tVt are suppressed by one order. Compared with branching fractions of the Be —> rjelv and 
Be J/ipW decays, the branching ratios of Be {he, Xc0,i,2)li^ are smaller by a factor of 2 — 10. The polarizations 
oi Be ^ helD and Be — » Xci''^ are dramatically different: it is very large for Be — > heW but very small for 
Be — > XcO^^- Most of our predictions on the semileptonic Be decays are comparable with results in the literature for 
the Be {he,XcO,i,2)l'^ decays. These results will be tested at the ongoing and forthcoming hadron colliders. 
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APPENDIX A: SOME SPECIFIC RULES IN THE p INTEGRATION 



When performing the p~ integration, we need to include the zero-mode contribution. This amounts to performing 
the i nteg ration in a proper way in this approach. To be more specific, for p[ under integration we use the following 
rules 



10 



111 



PluPluPl a. 



P'i^.p'iu = 9^..A^? + P^.P.A^2^ + {P^Sh + q^.Pu)A''p + g^g.Af , 



qa + Qt^aqv + guaqt,)A^2 
+P^,P,Pc.Af + (P^P.fc + P^g.P^ + q^^PMAf 

+ {qf,qi,Pa + qf.Pi^qa + Pf,qi,qa)A'i?^ + qf,q^qaA''Q\ 



9m 



4^'^2- 



P'luP'iv^'i = a^iuA'^^^ Z2 + qt,q^ 



AfZ2 + 2i:^A«Af)' 



(Al) 



where the symbol = reminds us that the above equations are true only after integration. In the above equation, A 
are functions of a;i,2, P±, p'± • q± and q^. Their explicit expressions have been studied in Ref. [3, EH- 



(i) 



N[ + mf -ml + {l^ 2xi)M'^ + {q^ + q ■ P) 



Pi_ ■ 9-L 



A[' 



A 



(2 



A 



(2 



A 



(3 



A 



(3 



Xi (1) _ (1) P^-q± 

T' ^2 -A, 



,2 iP'± ■ <1±? 

'Pi- 3 — ' 



- {A[ 



{A 



1 



tA 



(2) 



2 ) ^2-1 ' 



A 



(3) 



(l)^2 
1 J ' 

A[''Ar\ 



A 



(2) 



.(1) 4(1) 
^1 ^2 ' 



A 



(3) 



.(1) 4(2) 
^2 ^1 I 



.(1) .(2) 
^1 ^2 ' 

l(l) 4(2) 



A 



(3) 



.(1) 4(2) 
^2 ^2 ' 



(3) 



4(1) 4(2) 
^1 ^4 J 



4^-^^ A^^' H 4(1) 4(2) 

q 



(A2) 



We do not show the spurious contributions in Eq. (jA2p since they are numerically vanishing. 



APPENDIX B: EXPRESSIONS OF Be ^ A,T FORM FACTORS 



In this appendix, we collect the analytic expressions of Be he, Xi,2 form factors in the covariant light-front quark 
model. 

(Bl) 



V^^^^iq') 



1 



N, 



M' - M" 167r3 



dx2d pj^ 



X2N[N[' 



9 // 

q^w'{^ 



2xi{M' 



'2 ' Af^2) 



11 



f-q.P- 2{q^ + q ■ P)^^^ - 2{m[ + m'/)K - ^2) 



X2q^ +p'^ ■ q± 



xiq^w'lj^ 



p'_l-p"l + {x\m2 + X2m[){xim2 + X2m'[) 



2M 



/ dX2dy^^^^I^ 

J X2N{Ni' 



+^2(3^^ - 2Af ^ - 1) + ^ [^1(9' + 9 • P) - 2M'2 - 2p'^ • - 2m'i(m2 - m'/) 



2m2(m'i - m2)](^f ^ + - 1) + g • P 



(1) 



g2 + g4 



(44^) -3)1. 



(B2) 



(B3) 



(B4) 



(M' - M") 
1 



167r3 



M' - M " IGtt^ 



^X2N[N!l 

h' h" 

dx2dy^^x^ 

^X2N[N'{ 



+2x2m'^q ■ P + 2m2q^ - 2xiTO2(M'2 + M"^) + 2{m[ - m2){m[ - m'(f 



+S{m'i — 7712) 



,2 , {p'± ■ <i±y 



P± 



+ 2{m[-m'l){q'^ + q-P) 



p± ■ qi- 



= (M'-M")#^ f dx2d'p'^^^i^ 



{xi - X2){x2m[ + Xim2) - [2x\m2 - m'{ + {x2 - Xi)m'-i\ ^^ 



(B5) 



(B6) 



(B7) 



M' - M". 
2M" 



-V, 



M' + M" 



2M" 2 W I 2M" 167r3 



h' h" 

ax2d pj^- 



X2N[N[' 



2{2xi — 3)(a;2m'i + Xim2) — 8{m'i — m2) 



p'l , .{p'^-qi.? 



^ + 2 
r 



-[(14 - 12a;i)mi + 2m'/ - (8 - 12a;i)m2]^^^|. 



(B8) 



mv = 777 2 = m.c, in the case of the outgoing p-wave charmonia, so is zero, leading to the vanishing 
term of — ^ in the form factors. 



h^^^'^iq^) = fdx2dy^^^^{x2m[+x,m2 + {m[-m'{)?^ 

lOTT'^ J X2N[N'{ y r 

2 



,2 , (p1 • g^)' 
p^ "I 1:2 



+ 



K - 0(4'^ + + {< + m; - 2m2)(A^"^ + 



1(2) 



(2) , ^(2). 



mU4^^ + 4^ V ^7(24'^ + 24'^ - 4^^) 

HJrp 



(B9) 



167r3 



^a;2iV(iV{' 



2a;i(m2 - m\){M'^ + M^^) - Axim'iM'^ + 2x2m[q ■ P 



12 



+8(m'i - 1712) 



,1 , (Pj • 



+ 2(m'i+m'/)(9^ + g-P) 



Pi • 11. 



2xi{M'^ + M'a^)-q^ -q-P- 2{q^ + q ■ P) 



p± ■ q± 



—2(m'i — mi){m'i — 1712) 
Nr f , .0 , h'^h" 



dx2d pj_ 



2{A^^^ + A«)[m2((z' -N[- N'l - mf - mf ) 



-m[{M"^ - iV{' - mf - mj) ~ m'l{M'^ -N[- mf - ml) - 2m[m'lm,2] 

+2{m[ + m'/)(4''^2 + ^^i'^) + 16(m2 - m;)(Af ^ + A^^^) + 4(2m; - m'/ - ma)^^ ^ 



^ ( + M"2 - ^2 ^ 2(m; - m2)(m'/ + m2)](2Ap) + 2^'^ - A^^^) 



l(3) /l(2)^ 



(BIO) 



dx2cpp'±_ ^ T I (xi - X2){x2m\ + Xim2) - [2xim2 + m'( + {x2 - xi)m[]^^ — 

X2N[N[' [ q 



167r3 

X2q'^ +p'j_ ■ g_L 



9 // 

X2q^wj. 



P'l. ■ P" + (a;im2 + X2rni)(xim2 — X2m'{) 



+ 



167r3 



/ dx.dV^ |8(m2 - m;)(Af + 2Af ^ + 4^)) - 2m;(^« + A«) 

J X2N[N'{ y 



+4(2m'i - m'i - m2)(4^^ + 4'^) + 2(m'i + m'/)(4'^ + "^^f^ + 4'^) 

2[M'2 + - g2 + 2(m; - m2)K + m2)](4'^ + "^^f + 4'^ - 4'^ - 4'^) 



2 



1«) 



(3^) = - 



+[^2 _ ^/ _ _ + rn'/)'](^f + 2Af^ + 4^^ - A^^^ - a\ 

2{2x\ — 2>){x2rr}!y + X\m2) — 8(mi — 1x12) 



(Bll) 



/ dx2d^pj^- 



167r3 



"a;2iV(A^(' 



+ ^4 



-[(14 - 12a;i)m'i - 2mi - (8 - 12a;i)m2] ^-^ "/"^ 



( [M'2 + M"2 - g2 ^ 2(m; - m2)(m'/ + m2)](^r^ + 4^ - 4" 



(1)^ 



+^2(34'^ - 24'' - 1) + +q-P)- 2M'2 - 2pi • q^ - 2m; (m'/ + m,2) 



1(2) 



-2m2(m'i - m2)]{A^i^ + 4^^ -l) + q-P 



p'l , iPi-Q±y 

,2 



(44^^-3) 



dx2d^Pj_- ^ ^ 



q^ 

8(m2 - m'i)(4^^ + 24'^ + 4'^) - emU^l'^ + 4'^) 



167r3 J - '-X2N[N[ 
+4(2m'i - m'l - m2)(A^^^ + 4^') + 2(3m'i + m'/ - 2m2)(A<^^ + 24^' + 4^^) 



+- 



2[M'2 + M"2 - g2 + 2(m; - m2)(m'/ + m2)](4'^ + 24'' + 4'' - 4'' - 4'0 



1(3) 



1(3) 4(2) 



l(2)^ 



+2Z2(34'^ - 24'^ - 4'^) + 2^:^(64')^^) - 64') + l(Af ))2 - )) 
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+ - 2Af'2 + N[- N[' - (to'/ + ma)^ + 2{m[ - m^fX^A^^^^ + 2Af + - - (B12) 
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